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Abstract. We use the parabolic and hyperbolic equation with fractional time 
derivative to describe the subdiffusion in a system with thin membrane. We 
find the Green's function and solutions of the equation for the system where the 
homogeneous solution is separated by a thin membrane from the pure solvent. 
The solutions were obtained for two boundary conditions where the ratio of the 
concentrations at the membrane surfaces does not change in time and where the 
flux flowing through the membrane is proportional to the concentration difference 
between membrane surfaces. We discuss the difference between the solutions 
for parabolic and hyperbolic subdiffusion equations obtained for both boundary 
conditions. 
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1. Introduction 

The subdiffusion is usually denned as a process where the mean square displacement 
of a particle (Ax 2 ) is a power function of time 

where the subdiffusion parameter a is less than one (0 < a < 1), D a is the subdiffusion 
coefficient measured in the units m 2 /s a pQ. The case of a = 1 corresponds to the 
normal diffusion. The subdiffusion is related to an infinitely long average time that a 
random walker waits to make a finite jump. Then, its average displacement square, 
which is observed in a finite time interval, is dramatically suppressed. The subdiffusion 
occurs in systems with complex internal structure such as gels or porous media. To 
describe the subdiffusion the nondinear differential equation of natural order derived 
on the base of Tsallis formalism [2j [3] or the normal diffusion euqation with diffusion 
coefficient which is assumed to be a power function of time D(t) = Dat"^ 1 [4] were 
used. The Green's functions obtained for these equations fulfill the relation (p}, but 
the physical meaning of the assumptions providing to the equations is not always 
clear. For example, it is difficult to explain the decreasing of diffusion coefficient in 
time in a homogeneous system. The subdiffusion linear equation with fractional time 
derivative, derived on the base of Continuous Time Random Walk formalism [TJH], 
has not got such disadvantages. 

The subdiffusion in a membrane system was recently studied experimentally and 
theoretically. The motivation of the study is that the understanding of subdiffusion 
in a membrane system can help to model transport processes in systems so different 
as living cells and membrane microfilters (see for example [6]). The system with 
a membrane can also be used to measure subdiffusion parameters by comparing 
theoretical and empirical concentration profiles of substances of interest [7] . To model 
a transport process in such a system the parabolic subdiffusion equation was applied. 
However, the parabolic normal diffusion and parabolic subdiffusion equations give the 
solutions which posses the 'unphysical' property. Namely, for spatially unrestricted 
system the Green's function G(x,t;xo) (which is a probability density of finding a 
particle in the position x after time t under condition that at the initial moment 
the particle was located in xq) have non-zero values for any x and t > 0. This fact 
can be interpreted as an infinite speed propagation of some particles. To avoid this 
'unphysical' property Cattaneo proposed the hyperbolic normal diffusion equation 
based on the assumption that the diffusion flux is delayed in time r with respect to 
the concentration gradient [8] . The Green's function of the equation is equal to zero for 
finite x — xq, so the propagation velocity of the particles is finite. In phenomenological 
way the hyperbolic subdiffusion equation can be derived by involving the fractional 
time derivative into a flux or continuity equation. In the paper [9] there was noted 
that the hyperbolic subdiffusion equation can be derived in three different manners 
and the equations obtained are not equivalent to each other. 

The delaying effect of the flux with respect to the concentration gradient seems to 
be stronger in a membrane system than in a homogeneous one, since the flux can be 
involved into boundary conditions at the membrane. So, the delaying effect can appear 
not only in the equation but also in the boundary conditions. As far as we know, the 
hyperbolic subdiffusion equation has not been applied to describe the subdiffusion in 
a membrane system yet. In our paper we compare the solutions of the parabolic and 
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hyperbolic sub-diffusion equation in a homogeneous system and in a system with a 
thin membrane. 

The problem of choosing a transport model in a membrane system is more 
complicated since one of the boundary conditions at the membrane is not set 
unambiguously. Two boundary conditions which are not equivalent to each other 
were used. First of them demands the constant ratio of the concentrations at the 
membrane surfaces [101 [UJ , the second one assumes that the flux is proportional to 
the concentration difference between membrane surfaces [121 03] . The qualitatively 
difference between them is manifested in the long time limit because the concentrations 
calculated for second boundary condition goes to a continuous function at the 
membrane, unlike than for the first one. 

In our paper we find the solutions of the hyperbolic equations for a system with 
a thin membrane for two boundary conditions mentioned above and compare them 
with the ones obtained from parabolic equation. We consider the system where the 
thin membrane separates a homogeneous solution from a pure solvent (we add that 
such a system was often used in experimental studies QTJ 03 El)- In °ur study 
we assume that the system is one-dimensional where the diffusion or subdiffusion 
parameter as well as the parameter of membrane permeability do not depend on time 
and concentration, the first one is also independent of space variable. 

The paper is organized as follows. In section [21 we present the phenomenological 
derivation of the hyperbolic equation for normal diffusion. We show plots of the 
Green's functions obtained for the long times for the homogeneous system without 
a membrane. In section [3] we present the hyperbolic equation and the Green's 
functions for the subdiffusive system. The boundary conditions at a thin membrane 
are derived in section [U Solutions of the hyperbolic equation for the system where 
the homogeneous solution is separated by the thin membrane from pure solvent are 
presented in section [5] To illustrate our considerations the function obtained in 
sections [21 [3] and [5] are shown in several plots. Analyzing the plots we discuss the 
properties of the solutions in section [5l 



2. Normal diffusion equation 

2.1. Parabolic equation 

It is well known that the normal diffusion equation 
9C(x,t) _ n 9 2 C(x,t) 

at dx* () 

with normal diffusion coefficient D (measured in the units ro 2 /s), can be derived 
phcnomcnologically by combining the first Fick's law 

J6*)— B»gsa, ( 3) 

and the continuity equation 

dC(x,t) _ dJ(x,t) 

dt Ox ■ U 

The Green's function is defined as a solution of the equation for the initial condition 

G(x,t;xo) = 5(x - xq), (5) 
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and boundary conditions appropriate for considered system. When the system is not 
spatially restricted, there is 

G{— oo, t; xq) = G(oo, t; xo) = 0, (6) 

and the Green's function reads 

The function ([7]) is different from zero for any x and t > 0. Utilizing the 
probability interpretation of the Green's function one concludes that some particles 
are transported with the infinite speed propagation. 

2.2. Hyperbolic equation 

To ensure the finite velocity of the particle propagation one assumes that the flux is 
delayed with respect to the concentration gradient 

J(x,t + r) = -D^^, (8) 

where r is the delay time. Assuming that the parameter r is sufficiently small, the 
left hand side of equation ([5]) can be approximated by the first two terms of Taylor 
series with respect to r 

J{x,t)+T aJ % t) =-D ™g; t) . (9) 

Applying the operator d/dx to equation ([9]) and taking into considerations the 
continuity equation Q one gets the hyperbolic diffusion equation 

d 2 C(x,t) dC(x,t) „ d 2 C(x,t) 
T ^i^ + ^i— = D ^x^- (10) 
We add that equation [|lUp can be derived from differential-difference equations with 
continuous time and discrete space variable [16] . The process can be interpreted as a 
process with 'minimal' memory which extends to one time step more than in 'ordinary' 
diffusion process described by parabolic diffusion equation. 

Equation (|10p ensures the finite propagation velocity of the particles v = yjD/r. 
In the limit t — > we get parabolic diffusion equation with infinite v. To solve equation 
(1101) we must take two initial conditions. Let us assume that one of them is 



dC{x,t) 



= 0, (11) 

t=o 



dt 

what means that at an initial moment the concentration does not aim at its change 
and is effectively changed after time r since the particle flux is not generated before 
this time. The second boundary condition reads as C(x,0) = f(x). 

2.3. Green's function 

We obtain the Green's function for equation (|10P solving it by means of the 
Laplace L[f(t)] — f(s) — J °° f(t) exp(— st)dt and Fourier F[g(x)] = g(k) = 
I^oc 9^ x ) ^ x p{ikx)dx transforms method for the initial conditions |5]| (with xq = 0) 
and (fTTj) . After simple calculations we get the Green's function in terms of Lapalce 
and Fourier transforms 

G(M;0)= , 1+ 2 T °, 2 - (12) 

s + ts z + Dk z 
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The inverse Fourier transform of equation (|12p reads 

\x\y/s 



G(x,s;0) 



• exp 



(13) 



2v Ds V VD 

The hyperbolic equation was derived on the assumption that we omit the terms which 
iclude the parameter r k , k > 1, in the Taylor series of the flux (see equation (J9])). We 
find similar approximation of equation (|13[) . namely 

1 /„ TyG \x\s \ /" |l|v / s'' 



G(cc,s;0) 



1 



2VDs V 2 2VD 
The inverse Laplace transform of equation (fT"4|) is 

1 \x\t \ ( x 2 \ \x\t 

4Dt ) 



exp 



G(x,t; 0) 



1 



exp 



2VnDt V ' 4tV£>y 
where the function / is defined as [T7] 

fv,p(t;a) = L^ 1 [s^exp (- 



AD 



1/2,1/2 



D 

x\ 
~D 



(14) 
(15) 



)] 



for a, (3 > 0. This function can be expressed by the Fox function H and reads 



1 



1 

= ^z2 



ttW 



t 



k=0 



1 1 

(l + v)/P 1/(3 

/ a n k 
k\T(-kP - v) \¥, 



1 



(16) 



O 0,2 



100 
50 
10 




Figure 1. The plots of the normal diffusion Green's functions for different values 
of parameter r given in the legend, here t = 500, D = 10~ 3 . 

The plots of function (TT5)) are presented in figure Q] for different values of 
the parameter r. As we can see, only relatively large values of r make the 
noticeable difference between the Green's functions obtained for the parabolic equation 
(represented by the solutions for t — 0) and the hyperbolic one. 
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3. Subdiffusion equation 



3.1. Parabolic equation 

The hyperbolic subdiffusion equation can be derived by analogy with the derivation 
of the parabolic one. There are few ways to find the parabolic subdiffusion equation 
in phenomenological way. In the following we consider two of them which are natural 
generalization of the derivation of the parabolic normal diffusion equation. In the first 
one it is assumed that the subdiffusive flux reads as 

d 1 ^ dC(x,t) 



J(x,t) 



-D a - 



(17) 



'at 1 -* dx ' 

where d^ L /dt a denotes the Riemann-Liouville fractional time derivative defined as 
QH US] (here a > 0) 

{t-u) a - l f{u)du, (18) 



RL 

dt- 



1 



T(a) 



and 



d% L f(t) d» dg£"/(t) 



(19) 



dt a dt n dt a ~ 

where n is the lowest natural number that fulfills n > a. The Laplace transform of 
the Riemann-Liouville derivative is 

-k- 



L 



dt a 



JOL- 



7(*) 



fe=0 



dt' 



-k-l 



(20) 



t=o 



where n — 1 < a < n. Combining equation (jTTJ) with equation ^ one gets the 
parabolic subdiffusion equation [TJ [S] 

d 2 C(x,t) 



dC(x,t) _ 



(21) 



dt ""dt 1 -" dx 2 
In general, to solve the differential equation with the fractional Riemann-Liouville 
derivative by means of the Laplace transform method, one should fix the initial 
condition for time derivatives of the fractional negative order (see equation ([20]) ). 
what is beyond of physical interpretation. However, this remark does not concern the 
subdiffusion equation ([21]) since for a limited function there is (see Appendix) 

d^C{x,t) 



dt a - 



= 0, 



(22) 



when < a < 1. The Laplace and Fourier transforms of ([2"T]) is 

sC(k, s) - F[C(x, 0)] = -D a s 1 ~ a k 2 C{k, s). (23) 

For the Green's function with the initial condition ([5|) we get F[C(x,0)] = 1, what 
leads to the form of equation ([2"3"]l obtained from the Continuous Time Random Walk 
formalism [1]. 

Another scenario provides to the subdiffusion equation consists in changing the 
time derivative of natural order to the fractional Caputo one in the continuity equation 
([!]) according to the formula d/dt — > Od^/df* , where 9 is a parameter which is involved 
to achieve an appropriate physical units. The Caputo fractional derivative is defined 
by the relation []]|] 



d a c f(t) 
dt a 



1 



T(n — a) 



(* 



dt' 



;f(u)du, 



(24) 
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L 



dg/(*) 



dt a 



n-l 



= s a f(s) - £ 



fc=0 



t=0 



where n — 1 < a < n. Thus, we get 

fl 0gC(:E,f) _ gj^t) 

In the following we take 
D a 

Combining ([3]), (f26|) and (|27|) one gets the subdiffusion equation 

dgC(x,t) _ d 2 C(x,t) 
dt a a Ox 2 ' 

Equation (|28p is equivalent to equation (|2I[) since their Laplace and Fourier transforms 
are expressed by 



(25) 

(26) 
(27) 

(28) 



3.2. Hyperbolic equation 

The hyperbolic subdiffusion equation can be obtained by introducing the time 
derivative of fractional order to equations (J4J) or (J9j). As was noticed in the paper [9] , 
where the Riemman-Liouville fractional derivative only was taken into considerations, 
it can be done in three different manners. Unlike in [9], we involve Caputo fractional 
derivative into the continuity equation ([JJ . We assume that the flux is given as follow 



(29) 



J(xt + r)--D d ° M 

J(X,t + T)- U aQtl _ a ^ . 

Similarly to the previous case, let us approximate the left hand of equation (|29[) for 
t <C t by the first two terms of Taylor series with respect to r 



J(x,t) + T 



dJ(x,t) _ dC{x,t) 



at "at 1 -" dx 

From equation ([30]) and equation ((4]) we get the hyperbolic subdiffusion equation 



(30) 



_d 2 C(x,t) dC(x,t) 



= D a - 



l-a 
RL 



d 2 C{x,t) 



(31) 



dt 2 dt "dt 1 -* dx 2 ■ 

Let us note that from equations ([9]), (|26|) and |27|) we get the hyperbolic subdiffusion 
equation with Caputo fractional derivatives 



_ d 1 c +a C{x,t) | dgC(x,t) 



D a 



d 2 C(x,t) 



(32) 



9*!+° dt a " dx 2 ■ 

Equation (|32|) is fully equivalent to equation (fBTj) since the Laplace and Fourier 
transforms of the equations are the same. 
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3.3. Green's function 

As previous, we take the initial conditions {5J (for x = 0) and (TTTj) to solve equation 
n"j) . After calculations we get 

1 + TS 



G{k,s;0) = 



S + TS Z 



the inverse Fourier transform of equation ()33[) is 



G(x,s;0) 



■ exp 



„a/2 



(33) 



(34) 



The hyperbolic equation was derived under the assumption that we take into 
account linear terms in the Taylor series of the flux (see equation ©). Let us perform 
similar approximation for equation (|f 3[) . what gives 

tW 2 \ X \s a \ ( \x\s a / 2 



G{x,s;0) 



1 

2jDls 1 - a / 2 ~ 



1 



|a;|s Q 
2JK 



exp 



(35) 



The inverse Laplace transform of equation (|35|) is 



G(x,t; 0) 



1 



fa/2-1, a/2 I t 



.C 



2^ 



/a, a/2 I 



I:r| 



/2,a/2 



(36) 




Figure 2. Hyperbolic subdiflusion. The plots of the Green's functions for 
different values of r, here t = 500, D a = 10 — 3 , a = 0.5. 

The plots of the Green's functions ([56]) are presented in figures Contrary to 
the normal diffusion case, the effect of delaying is hardly observed in the considered 



cases. 
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Figure 4. The Green's functions for t = 500, 1000, 1500, 2000, here D a = 10 -3 , 
a = 0.8, the dashed lines correspond to r = 100, continuous ones represent the 
Green's functions with r = 0. 



Hyperbolic versus parabolic equation with fractional derivative 



10 



4. Boundary conditions at thin membrane 

We denote the concentration and flux in the region x < x m as G\ and J\, and in 
the region x > x m as C2 and J2, respectively, where x m is the membrane position. 
Since the equation is of the second order with respect to the space variable, we need 
two boundary conditions in each of the region. Two boundary conditions demand 
finiteness of the solutions at x — > — 00 and x — > 00, two other ones are fixed at the 
membrane. The first of them is rather obvious and it assumes the continuity of the 
flux at the membrane 

Ji(x ml t) = J 2 (x+,i) = J(x m ,t). (37) 

However, the problem of fixing the second boundary condition at the thin membrane is 
not unambiguously solved. The missing boundary condition at the membrane has been 
chosen in two ways. The first one demands the constant ratio of the concentrations 
at two opposite sides of the membrane [7J [THl HJ 

C 2 (a;+,i)/Ci(x m ,i) = 7 = const. (38) 

In the second one the flux flowing through the membrane is proportional to the 
difference of the concentrations at the opposite sides of the membrane [Jj [T2l [13] 

J(x m ,t) = X [C 2 (x~ , t) - Ci (x+ , t)} , (39) 

where A is the membrane permeability coefficient. Below we consider the possibility of 
application of the boundary conditions ([38]) or (f39|) for the case of the system described 
by the hyperbolic normal diffusion or hyperbolic subdiffusion equation. 



4-1. Constant concentration ratio at the membrane 

In the Smoluchowski's papers (see for example [20) it was derived the boundary 
condition at the fully reflecting wall. Let the wall is placed at x m and the system 
occupies the interval (— oo,a: m ). Smoluchowski's original approach utilized the 
assumptions that the amount of the substance in the system does not change in time 

^- t J m C{x,t)dx = Q, (40) 

and the flux vanishes at —00. Integrating the continuity equation (|H) over the interval 
( — 00, x m ) and using the above assumptions one gets 

J(x m ,t) = 0. (41) 

Since equation ^ works in hyperbolic subdiffusion equation case, we choose the 
boundary condition (|41 [) at the fully reflecting wall for the system described by the 
hyperbolic equation. 

Chandrasekhar used the method of images to derive the Green's function for 
this system [21] . The Green's function can be interpreted as a concentration of large 
number of particles N (divided by N); the particles are located at point Xo at the 
initial moment t — 0. So, the Green's function can be treated as an instantaneous 
particle source (IPS) normalized to 1. Within the method one replaces the wall by 
additional IPS in such a manner that the concentration behaves exactly as in the 
system with the wall. Vanishing of the flux at the reflecting wall is achieved when one 
replaces the wall by the IPS located symmetrically to the initial point xq with respect 
to the wall. Then the instantaneous particle sources create fluxes of particles flowing 
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in opposite directions, which reduce each other at the point x m . Thus, one finds (for 
x < x m and x < x m ) 

G(x, t; x ) = G (x, t; x ) + G (x, t; 2x m - x ), (42) 

where Go denotes the Green's function for a homogeneous system (without the wall). 
Let us note that the Green's function equation (jl2")) leads to equation in any 
system where the flux fulfills the following relation 

'-§■ 

In the papers [TU] the method of images was generalized to the system with a 
partially permeable wall. Since the Green's function (|42|) works in the system 
with fully reflecting wall, where a transport is described by hyperbolic diffusion or 
subdiffusion equation, similar generalization can be performed in such a system with 
thin membrane. 

In the system with thin membrane the particles can pass trough the membrane in 
both directions many times. Let us assume that the membrane is symmetric and the 
probabilities of passing through it do not depend on the direction of particle's motion. 
To take into account selective properties of the membrane, we 'weaken' the additional 
IPS located at 2x m — Xq by factor a, which gives 

G(x, t; x ) = G (x, t; x ) + aG a {x, t; 2x m - x ). (44) 

Assuming that the flux is continuous at the membrane, we get for x > x m and xo < x m 

G{x,t;x ) = (l-o-)G {x,t;xo). (45) 

The functions (jSJ) and (gSJ) fulfill the boundary condition (J3SJ) with A = (1 - 
<t)/(1 + er). To interpret the parameter a let us note that the probability of finding 
a particle (starting from xq, where xq < x m ) in the region x > x m is equal to 
Pa = (1 — a) J Go(x,t;xo)dx for the membrane system, whereas the probability 
of finding the particle in this region for system without the membrane is equal to 
P<3 = J x Gq(x, t; xo)dx. Comparing the above equations we obtain a = 1 — P a j Po, so 
the parameter a can be interpreted as a probability of finding the particle in the 
region x < x m under condition that in the system with removed membrane the 
particle will be in the region x > x m . In other words, a is a conditional probability 
of stopping the particle by the membrane in unit time under condition, that in the 
similar system with no membrane this particle pass the position x m . Thus, a is the 
parameter controlling the reflection of particles by the membrane, the parameter 1 — a 
is the parameter of membrane permeability. The boundary condition (|38[) has simple 
physical interpretation: if N diffusing particles are going to pass through the wall in 
unit time, then o~N of them will be stopped by the wall whereas (1 — o~)N pass through, 
where a = (1 - A)/(l + A). 



4-2. Radiation boundary condition 

For the parabolic normal diffusion equation the radiation boundary condition (|39p was 
derived from the model with discrete space variable [12j as well as for the considerations 
performed in a phase space where the diffusion is described by the Klein-Kramers 
equation [13]. Equation (|39|) can be interpreted as the natural continuation of the 
Fick equation applied to the membrane. According to equation (|8|), we can generalize 
equation (|39|) as follows 

J(x m ,t + T) = \[C 1 (x-,t)-C 2 (x+,t)]. (46) 
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In the following we will see that the boundary conditions and (|46[) are not 

equivalent to each other. 



5. Solutions of hyperbolic subdiffusion equation for a membrane system 

0. We choose the initial 



C(x,0) 



(47) 



Let us assume that the thin membrane is located at x n 
condition as 

Co, x < 
0, x>0. 

The boundary conditions demand finiteness of the solutions at infinity 

C 1 (-oo ) t) = C , C 2 (oo,f) =0. (48) 

The first boundary condition at the membrane (|37|) ensures that the flux is continuous 
at the one, the second boundary condition at the membrane we take in general form 

hCi(0-,t) +b 2 C 2 {0 + ,t) + b 3 J(Q,t + T) = 0. (49) 

We note that the Laplace transform of the flux reads 

° X - Q dC(x,s) 



J(x, s) 



-D n 



(50) 



1 + ts dx 

The Laplace transforms of solutions for boundary conditions (|37|). (|48j> . (|49l and initial 
condition (|4"7) are as follows 



Ci(x,s) = 



C 2 (x,s) 



Cn 



Cn 



1 



h-b 2 - b?,^D^s 1 - a / 2 /y/T+T~s~ 



exp x 



bi 



: exp — x 




(1 + rs)s c 



(1 + Ts)s a 



(51) 



(52) 



s bx-b 2 - b 3 VD^s 1 ~ a / 2 /VTTTs~ ~' r \ V 
Below we find the solutions for two boundary conditions described in section [U 

5.1. Constant ratio of the solutions at the membrane 

Putting b x > 0, b 2 < and b 3 = in |g5J) we get equation ([35]) with 7 = -b 2 /h. The 
solutions are as follows 



Ci(x,t) =C 



1 - c/-l,a/2 I 



xt / -ac 

0'^^=/q/2,q/2 I * 



2 V / £'a' 



Do 



C 2 (x,t) =C o- 



/-l,a/2 I ^ 



a/2, a/2 I ^ ' 



.(54) 



where a = 1/(1 + 7). The plots of functions ([55]) and ([54")) are presented in figure [5] 

As we can see, the differences between the solutions obtained for the parabolic 
subdiffusion equation are very close to the solution of hyperbolic equation (even for 
the largest time r = 100). 
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Figure 5. The solutions calculated for the boundary condition H38H with 7 = 1.5, 
a = 0.9, D a — 5 X 1CT 4 for t = 500, 1000, 1500, 2000. Vertical lines represent the 
membrane, dashed lines correspond to r = 100, continuous ones correspond to 
t = 0. 



5.2. Radiation boundary condition 



Here b\ = —62 > 0, 63 < 0. The boundary condition takes the form of equation (|46[) 
where A = —61/63. To obtain the inverse transforms of equations (|51[) and (j52"|) we 
assume that ts C 1 (what corresponds to f > l/ r ) anc ^ we extend the transforms 
into the power series with respect to the parameter s. Achieving only the linear terms 
with respect to r we get 

*y 00 / rr^ — \ k 
<^0 



Cl (M)=Cb-^E(-^ 



2 ^ \ 2A 

At 



/fc(l-a/2)-l,a/2 I * 



2~/fc(l-a/2),a/2 ^5 ~ 7g= 



XT I —X 

^-y=/fe(l-Q/2)+a/2,a/2 I *5 



65) 



C 2 (M) 



Co 



E 

fc=0 



/Do 
2A 



fc r 



/fe(l-a/2)-l,a/2 I * 



fcr ( x 

+ "7r/fe(l-a/2),a/2 I * 



/ a: 

-Jk{l-u/2)+a/2,a/2 \ l > 



The plots of functions ([55| and lj5l)|) are presented in figures [5] and [7] 



Da 

(56) 



6. Final remarks 



We present here the solutions of the parabolic and hyperbolic subdiffusion equations 
for the homogeneous system and for the membrane one. The solutions were found 
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under assumption that we take into account the terms linear with respect to the 
parameter r. We applied two different boundary conditions at the membrane. Our 
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considerations are illustrated by few plots presenting the solutions for both of the 
boundary conditions. The plots were prepared for the parameters which values are 
of the order of the ones already found for real systems on the basis of experimental 
results [7]. The detailed remarks extracted from the plots are not fully conclusive, but 
it suggested few regularities, which - in our opinion - are general. They are as follows. 

(i) For the boundary condition (|38|) 

• The solutions at the membrane do not change in time and read as C\ (0~ , t) = 
(7Co)/(l+7), C2(0 + ,£) = Co/(l+7). This property seems to be 'unphysical', 
but let us note that the solution obtained for the system without membrane 
(for which 6 X = — b 2 and 63 = 0) with initial condition (|¥7)l are constant for 
x = and reads as C(0,<) = Co/2. 

• The delay effect does not occur at the membrane, consequently it is weak at 
the membrane neighborhood. 

(ii) For radiation boundary condition (|46|) 

• The concentration difference between the surfaces decreases in time (see 
figure [6]). From equations (fSTj) and (|52|l it is easy to see that Ci(0, t) — > 
6*2(0, t) when t — > 00, since the long time limit corresponds to the limit of 
small s. 

• For A ~ 10 _1 the membrane loses its selectivity. 

(iii) The main qualitative difference between the above boundary conditions is 
noticeable in the long time limit as boundary condition (|39[) leads to the solutions 
of hyperbolic subdiffusion equation which are going to the continuous function at 
the barrier, unlike than the solutions obtained for ([38]) . 

(iv) In all cases the delayed effect is connected with the subdiffusion parameter a (this 
property is clearly seen in figure Q]- figure [3] for the Green's functions). When a 
increases, the delaying effect is stronger. 

(v) For large times the delaying effect is negligibly small. Let us note that for large 
time the term r/t vanishes (it corresponds to the limit ts — > in equations (f5Tj) 
and (j52"l) ). 

Here the question arises: why the subdiffusion hyperbolic equation has not been 
applied to describe the experimental results in the subdiffusive membrane system, 
despite of proper 'physical quality' of the equation? Analyzing the plots [2] - [7] 
we conclude that in considered cases there is no reason to apply the hyperbolic 
subdiffusion equation instead of the parabolic one. The difference between the 
solutions is so small that both of them would certainly be laid within the error bars 
of the experimental concentration profiles. The order of values of the subdiffusion 
coefficient D a taken into calculations agrees with the ones obtained experimentally 
for sugars in agarose gels [7] if as unit of time 1 sec is chosen and 1 mm is the unit 
of space variable. In these units the value r = 100 is certainly too large, nevertheless 
these differences are rather hard to observe, for smaller values of r these differences 
are smaller. 

There is a problem with choosing the boundary condition at the membrane. 
Seemingly there is no problem with choosing the condition since the real system is 
limited by external walls and the concentration goes to equilibrium functions which is 
continuous at the membrane. This property posses the radiation boundary condition 
only. However, experimental study performed on two membrane system show that 
the concentration profiles have the scaling property, which is limited to the theoretical 
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solutions obtained for the boundary conditions only [T5J [55] . Namely, changing 
the variables according to the relation (x,t) — > (p a x,pt), p > the experimental 
profiles does not display noticeable changes. So, if the particles flowing through the 
membrane do not 'feel' the presence of external walls of the system, the boundary 
conditions (|38|) can be used. 

Although our conclusion is rather odd in respect to the membrane system, there 
are systems where the solutions of the hyperbolic and parabolic equations considerably 
differ from each other. Such a situation occurs for the boundary conditions where the 
concentration of the particles oscillates with high frequency, as for example in the 
problem of impedance spectroscopy [23J . 
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Appendix 

Here we proof the relation ([55)1 . Let us assume that \C(x,u)\ < A for u 6 (0,i\ and 
C(x, t) ->• when t -> 0. Then, 



From the above equation we get equation ([551) . 
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